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A Prototype Novel Sensor for Autonomous, Space Based Robots 


I. Introduction 

The goal of this program was to develop new sensing capabilities for 
autonomous robots operating in space. Information gained by the robot using these 
new capabilities would be combined with other information gained through more 
traditional capabilities, such as video, to help the robot characterize its environment 
as well as to identify known or unknown objects that it encounters. Several sensing 
capabilities using nuclear radiation detectors and backscatter technology were 
investigated. 

The result of this research has been the construction and delivery to NASA of a 
prototype system with three capabilities for use by autonomous robots. The primary 
capability was the use of beta particle backscatter measurements to determine the 
average atomic number (Z) of an object. This give the robot a powerful tool to 
differentiate objects which may look the same, such as objects made out of different 
plastics or other light weight materials. 

In addition, the same nuclear sensor used in the backscatter measurement can 
be used as a nuclear spectrometer to identify sources of nuclear radiation that may 
be encountered by the robot, such as nuclear powered satellites. A complete nuclear 
analysis system is included in the software and hardware of the prototype system built 
in Phase II of this effort. 

Finally, a method to estimate the radiation dose in the environment of the 
robot has been included as a third capability. Again the same nuclear sensor is used 
in a different operating mode and with different analysis software. Each of these 
capabilities are described below. 


II. The Prototype System 

The system consists of a custom made, highly sensitive CdTe nuclear sensor, a 
small isotopic source of beta particle radiation, a special collimator/detector package, 
a preamplifier, a personal computer, a multichannel analyzer, needed interfacing 
electronic circuitry, and all of the needed software. 

I 

A. Description of System Hardware 

1. CdTe sensors 

CdTe is a semiconductor material which can be used to fabricate very sensitive 
nuclear spectrometers which can operate at and above room temperature. CdTe 
was chosen because it high sensitivity allowed the entire sensing head to be small and 
because, by using CdTe no cryogenic coolant was necessary. In addition, these 
sensors have been used on previous space missions and have demonstrated there 
long term "space-worthyness". 

2. The Beta Particle Source 

During the program, a method was developed that allowed for the accurate 
determination of the average atomic number (Z) of a target independent of the 
density or chemical makeup of the source. This determination was done using beta 
particle backscatter. Beat particle scatter off electrons in a target. The number of 
electrons in the target is directly related to the atomic number of the constituent 
atoms and of the density. The amount of particles scattered is related to the density 
and the spectrum of the backscattered particles is related only to the average Z. 
Thus, by measuring the spectrum, the average Z can be determined. In principal, the 
density can be determined also by measuring the flux of backscattered particles. 
However, this is very difficult to calibrate and is a strong function of the distance 
from the source, the geometry of the measurement and to some extent the size and 
shape of the target. The backscatter spectrum is not affected by these factors. 

Strontium-90 , ^r, was chosen as the beta particle source because its long half 
life allows for stability and low maintenance and its high maximum energy, 2.7 MeV 
permits deeper interrogation into targets than would be possible with lower energy 
sources. 

Figure I shows the results obtained. Clearly, plastics and elements are 
distinguishable by there atomic number. The relationships held for compounds also. 




3. The sensor head 

Figure 2 shows a schematic of the sensor head that was developed for the 
prototype. Figure 3 and 4 show photographs of the unit. The primary design 
considerations were that it must 1) be compact and light weight, 2) be safe to 
operators in the vicinity, and 3) permit accurate measurement of the target of 
interest. The design satisfies all of these requirements. - 

The head consists of a Delrin plastic sensor head, a CdTe detector, a ^Sr 
isotopic source contained in a graphite holder, and a lead shield located directly 
between the source and detector. This head attaches directly to a compact PCT-3 
preamplifier manufactured by RMD, Inc. The preamplifier plugs into the interface 
ports on the personal computer. 

D. Description Tiie Computer System and System Operation 

The prototype makes use of a standard, IBM-type, portable computer. A 
commercial multichannel analyzer, "The Norland 5000" board by Norland Corp., has 
been added. All of the user interface software has been written in BASIC. Three 
capabilities are accesses from the keyboard through the menu screen that appears on 
boot up. 

1. Measurement of Z Effective 

This primary capability provide the robot with the average atomic number of 
the target. With the isotopic source provided, the system will measure the spectrum 
backscattered from the target and determine Z eff . Z cff is close to, not exactly the 
simple arithmetic average of the constituent elements but related to the atomic 
numbers of the them by the following equation: 

-7 =I(^/AJZ- Mi)Zl 

eur 


where Wj is the weight fraction of the i 1 ^ 1 element, Aj is its atomic weight and Zj is its 
atomic number. 

This function of the prototype is selected from the keyboard by first choosing 
the Zeff function from menu on the screen, as described in the attached operators 
manual. 

2. Dose Rate Measurement 

The gamma ray radiation dose rate dose at the location of the robot is 




Fig. 3 • Photograph of the prototype backacatter sensor with the front 
faceplate missing. Circular object is the front of the beta sensitive 
CdTe detector, oblong object is the lead x-ray shield, smaller hole is 
for the beta source. 



Fig. y . Photograph of the adjustable and fixed sensor heads with the 
faceplate collimators in place. 
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determined using an empirical method based on the response on the CdTe detector 
to gamma rays. The detector measures the average energy of many rays striking it, 
the average energy deposited in the detector is measured and the "tissue equivalent" 
dose is estimated using a predetermined calibration curve. Experiments performed 
with this detector and software showed the dose rate could be determined with about 
a 20% accuracy over the energy range of 100 keV to about 3 Mev. 

Measuring the dose rate is achieved by selecting "Personal Dosimeter" at the 
menu screen. 

3. Nuclear Spectrum Analyzer. 

The third capability included in the prototype is nuclear pulse height analysis. 
Incorporated into the prototype is a complete Norland 5000 multichannel analyzer. 
This feature is accessed by choosing "Norland" from the menu. The details for 
operating this analyzer are described in detail in the operators manual which 
accompanies the prototype. 

III. Summary 

This SBIR program has been successfully completed and has met the technical 
goals of developing new sensing capabilities for autonomous space-based robots. A 
complete and fully functional prototype has been assembled, tested and delivered to 
NASA at JPL in Pasadena. This prototype demonstrate several possible capabilities 
which should provide useful senses to robot designers for future generations of 
robots. 
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I. INTRODUCTION 


numhJr'ry ' ? T*” 11 '" inslrumem “ designed lo provide an accurate measurement of effective atomic 
number (Z clf ) ot a target in under one minute at a distance of 4 cm. from the sensor head. Z cff is defined 


Zerr - S(wi/A,)Zj 2 /2(wi/Ai)Zj 


atomi^number ofthe^f 1 * EST* * elemeM> * = ° f ,he “ «•««« * = «* 

Iilb pr °' 0type ™ nsists 0|: 0) a Delrin pUstic sensor head, (2) CdTe beta detector, (3) 
Sr/ Y beta source m a threaded graphite source holder, (4) standard PCT-3 preamplifier (5) portabl l 

andX tl sto°r mPUler | ^ f* mdU , deS SpCCially desi 8 ned expansion boards for providing detector bias^oltage 
* nd ,or itor, . n S ™ d dialyzing the data, (6) cables. The entire system is compact, powered by 110VAC and 
very convenient to operate in the laboratory environment. P by liUVAL and 

II. INITIAL SET-UP 

Since the computer has been provided with the expansion boards already installed the onlv 
interconnection to be made is to connect the preamplifier cord to the 9-pin connector^ m the rw of the 

runs between " AM P O U a a d d j 1 > u n . i t is worth checking to see that the short BNC cable which 

after shipping. 1 ^ Und ” MCA 1N " ° n the sec °" d board is still firmly connected 

the n'Wn'anifl^r T '"‘T™’ Sett ' ngS ° n the am P Iifier expansion board which can be used to adjust 
the g.nn and detector bias, these settings have been optimized at RMD and, in general, should not be 

III. OPERATION 

lisced below" "* pr ° gram ' ,ypc " ZEFF - The raain menu allows lh <= •»« •» select from ,hc five options 

[ ] INSTRUMENT SETTINGS 

[ ] BACKGROUND ROUTINES 

[ ] CALIBRATION ROUTINES 

[] MEASUREMENT ROUTINE 

[] EXIT THE PROG RAM 

A. Instrumernt Sirt-riNcs 

and c^lihr.n Pti ° n 3,1 — ° nC ,0 .* CI the Spectrum Re fi'°n Of Interest (ROI), background, data collection, 
setting.! acc l ulsltlon periods. Normally the ROI settings should not be changed from the defauh 

MENU° Ce thC INSTRUMENT SETTINGS have been chosen the operator should return to the main 

B. Background Routines 

Before using either the calibration or measurement routines a background spectrum must be stored 
Therefore the BACKGROUND ROUTINES should be chosen from .he main MENU. 


«» i— - in .he sensor heeO. 

C . Calibration Routines 

Values tor calibration constants A R .lnr-i r* • 

target. Therefore one of the options from the cal.W 10 ° rder to determine Z* of a given 

the MEASUREMENT ROUTINES. ‘ ,0n m ° de must be exercised before proceeding to 

Choosing the option "USE DEFAULT CAI TRR attam" . . „ 

which were preset at RMD. Dependine on the mrf i N • automat,ca,| y sets A . B . and C to values 
may nul yield accurate results P |„ ori f er to J'' Jr “Pmmienml arrangement, these values may or 

"RECALIBRATE " option should be c 1° T a new experimental set-up, the 
tire required for calibration. For each calibration 'Tm ° f l hr “ lareels wllh knowrl a,omic numbers 

and then a spectrum taken. After h back cm, 1? ‘"T" ^ nUmber Sh0uld bc ‘"tered 

measured, calculated values for A B mid C are dh„| Pe ?T„ ' aS ' Calibra,i °" < a ^' been 

using the determined values lor A,' b' and C ar P 7 ,! E ‘ eC " Ve <"° raic num ber values, calculated 
calibration. When the calibration is complex the “Z ‘° CValUi “ e ‘ hC aCcura<: >' of ,he 

Sy ' ,e The"°T-, ready '° de ‘ Crmine elfec,ive ’““"lie -'umber of un^Zn laZf “ “ "* mai " ™ enU ' ^ 
rccalibration ^ through the 

Using this option allows one to directly inZh^eTfo™ CONSTANTS" option should be chosen. 

numbers of ^ ^ ,he of atomic 

calibration samples with Z u , ^ be — 6 a "“ 3 °. •*«" 

appropriate. ' y th y |ene )» 13 (Aluminum), and 29 (Copper) would be 

D. Measurkment Routine 

spectra and ca!ibrat^^^hLe^e™kw Uno" should . be chosen aft er appropriate background 
software displays TOTAL COUNTS IN ROT a vcdaSc'? 1 otac 9 u "'ing data and calculating Zc frt the 
ATOMIC NUMBER, CHANNEL NUMBE ^. a nd EFFECTIVE 

updated at the completion of the calculation. * pertormed * these numbers are cleared and then 

E. Exit Tiil: Program 

the computer is tmSeTte t^ Za Sd ‘° D0S ‘ Wben in D0S ' 
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■r ************************** 

' * ********************** 

* 

' * ZEFF.BAS * 

' * 

* 

' * ( RMD Inc. ) * 

• * 

* 

• * 


' *** 




DECLARE FUNCTION chanread& (A%, B%, C%) 
DECLARE SUB erasepha (A%, B% , C% , D%) 
DECLARE SUB setlld (lld%) 

DECLARE SUB setuld (uld%) 

DECLARE SUB setoffset (offset%, gain%) 
DECLARE SUB acquire (A%, B%, C%) 

COMMON SHARED TIME% 


' READ DATA A CHANNEL AT A TIME 
7 ERASES PHA MEMORY 
7 SETS LOWER LEVEL DISCRIMINATOR 
7 SETS UPPER LEVEL DISCRIMINATOR 
7 SETS DIGITAL OFFSET AND CONVERSION GAI 
7 CONTROLS ACQUISITION TIME 

7 ACQUISITION TIME 


DIM backgroundarray& (1024) 7 
DIM dataarray& ( 1024 ) ' 
DIM processedarray ! (1024) ' 
DIM Z (20) / 
DIM E (20) / 

DIM instdata(20) / 


ARRAY WHERE BACKGROUND IS STORED 
ARRAY WHERE DATA IS STORED 
ARRAY FOR PROCESSED DATA 
CALIBRATION Zeff ARRAY 
CALIBRATION <CH#> ARRAY 

INSTRUMENT SETTINGS 


DATA = 0 
BKGND = 
CAL = 0 
KEY ( 12 ) 


instdata 

instdata 

instdata 

instdata 

instdata 


0 


ON 


(10) 

= 15 

(12) 

= 800 

(14) 

= 100 

(16) 

= 15 

(18) 

= 30 


7 FLAG, 0 = NO DATA ACQUIRED 
7 FLAG, 0 = NO BACKGROUND INFORMATION 
7 FLAG, 0 = NO CALIBRATION INFORMATION 


7 ROI LOW 
7 ROI HIGH 

7 BKGND MEASUREMENT PERIOD (TYME) 
7 DATA MEASUREMENT PERIOD (TIME) 

7 CALIBRATION MEASUREMENT PERIOD 


A = -.2060272 
B = 14.07163 
C = 27 . 1889 


7 DEFAULT VALUE FOR QUADRATIC COEFF. 

' II 


A 

B 

C 


' * * 
BOX: 


* * * * * 


' ** 
/ 
/ 


*************** 

*************** 

* 

*** CONSTRUCTS 
* 


ON SCREEN DISPLAY BORDERS ************** 


* 

* 

* 


************** 

************** 

************** 


FOR V = 1 TO 23 

IF V = 1 OR V = 6 OR V = 21 OR V = 23 THEN 
FOR H = 1 TO 80 

LOCATE V, H: PRINT "e" ' HORIZONTAL LINES 


END IF 


NEXT H 
FOR H = 


IF V 
IF V 
IF V 
IF V 


1 TO 80 STEP 79 

LOCATE V, H: PRINT "n" ' VERTICAL LINES 

NEXT H 

1 THEN LOCATE 1, 1: PRINT "e": LOCATE 1, 80: PRINT "£" 

6 THEN LOCATE 6, 1: PRINT "a": LOCATE 6, 80: PRINT " 

21 THEN LOCATE 21, 1: PRINT "a": LOCATE 21, 80: PRINT 
23 THEN LOCATE 23, 1: PRINT "a": LOCATE 23, 80: PRINT " 

NEXT V 

LOCATE 2, 38: PRINT "ZEFF" 

LOCATE 3, 20: PRINT "TO MEASURE EFFECTIVE Z BY BETA BACKSCATTER" 

LOCATE 5, 3: PRINT DATE$ ? »" 

LOCATE 5, 71, 0: PRINT TIME$ ; "" 

*************************** ^^^^^^^^^^^^^^^^^ 

ACKNOWLEDGEMENT: ' ************************************************************* 

' * 

' *** PRINTS ON-SCREEN ACKNOWLEDGEMENT *************** 
' * 

9 **************************************************** 


LOCATE 

8, 

23, 

0: PRINT 

LOCATE 

12, 

24 

PRINT " 

LOCATE 

14, 

34 

PRINT " 

LOCATE 

16, 

26 

PRINT " 

LOCATE 

18, 

35 

PRINT " 

LOCATE 

20, 

33 

PRINT " 

LOCATE 

22, 

29 

PRINT "] 

DO 

LOCATE 5, 71, 


617-926-1167" 


PRINT ; TIME$; "" ' update real time 

— — ^ * "** j. j. ii i\ij i y — ’ ' • • 

GOSUB FLUSH 


10, 

11: 

PRINT 

12, 

11: 

PRINT 

14, 

11: 

PRINT 

16, 

11: 

PRINT 

18, 

11: 

PRINT 


[ ] INSTRUMENT SETTINGS" 
[ ] BACKGROUND ROUTINES" 
[ ] CALIBRATION ROUTINE" 
[ ] MEASUREMENT ROUTINE" 
[ ] EXIT THE PROGRAM" 


CO 


LOCATE 22, 14: PRINT " 


USE <UP> OR <DOWN> ARROWS THEN <ENTER> TO SELECT." 


D1 = 18 

MLIM : ' *** MANAGES THE LIMITS OF TRAVEL OF THE CURSOR ************************* 

IF D1 > 18 THEN D1 = 10 ELSE IF D1 < 10 THEN D1 = 18 
LOCATE Dl, 13, 1, 0, 13 

D1 = CSRLIN 

A$ = INKEYS: IF A$ = CHR$(13) THEN GOTO MCURSVAL ELSE IF A$ = "8" THEN D 
IF A$ = "2" THEN Dl = Dl + 2 : GOTO MLIM 
IF LEN ( A$ ) = 2 THEN A$ = RIGHTS (A$, 1) ELSE GOTO MLIM 
IF A$ - "H" THEN Dl = Dl - 2 : GOTO MLIM 
IF A$ * "P" THEN Dl - Dl + 2: GOTO MLIM 
GOTO MLIM 

Dl = Dl + 2: GOTO MLIM 

MCURSVAL: ' *** FINDS THE VALUE OF THE LINE AT THE 

Dl - CSRLIN 

IF CSRLIN = 10 THEN GOSUB INSTRUMENT 
IF CSRLIN - 12 THEN GOSUB BACKGROUND 
IF CSRLIN =* 14 THEN GOSUB CALIB 
IF CSRLIN = 16 THEN GOSUB MEASURE 
IF CSRLIN = 18 THEN BEEP: CLS : END 
GOTO MENU 


CURSOR LOCATION ************* 

' INSTRUMENT SETTINGS 
' BACKGROUND ROUTINES 
' CALIBRATION ROUTINE 
' MEASUREMENT MODE 
' EXIT FROM THE PROGRAM 




' * 

' *** INSTRUMENT SETTINGS MENU ******************************* 
• * 

GOSUB FLUSH 


LOCATE 7, 22: PRINT "*** 
LOCATE 10, 3: PRINT " ( ] 

LOCATE 12, 3: PRINT " [ ] 

LOCATE 14, 3: PRINT " [ ] 

LOCATE 16, 3: PRINT " [ ] 

LOCATE 18, 3: PRINT " [ ] 

LOCATE 20, 3: PRINT " [ ] 

LOCATE 22, 14: PRINT " 

Dl = 20 


INSTRUMENT SETTINGS ROUTINE ***" 

SET ROI LOWEST CHANNEL"; TAB(50); instdata(lO) 

SET ROI HIGHEST CHANNEL"; TAB(50); instdata(12) 

SET BACKGROUND ACQUISITION PERIOD"; TAB (50); instdata 
SET DATA ACQUISITION PERIOD"; TAB(50); instdata ( 16) ; 
SET CALIBRATION ACQUISITION PERIOD"; TAB(50); instdat 
RETURN TO MENU " 

USE <UP> OR <DOWN> ARROWS THEN <ENTER> TO SELECT." 

' SETS CURSOR TO PROPER 


ILIM: ' *** MANAGES THE LIMITS OF TRAVEL OF THE CURSOR ************************* 

IF Dl > 20 THEN Dl = 10 ELSE IF Dl < 10 THEN Dl = 20 
LOCATE Dl, 5, 1, 0, 13 

Dl = CSRLIN 

A$ = INKEY$ : IF A$ = CHR$(13) THEN GOTO ICURSVAL ELSE IF A$ = "8" THEN D 

IF A$ = "2" THEN Dl = Dl + 2 : GOTO ILIM 

IF LEN (A$) = 2 THEN A$ = RIGHT$(A$, 1) ELSE GOTO ILIM 

IF A$ = "H" THEN Dl = Dl - 2 : GOTO ILIM 

IF A$ = "P" THEN Dl = Dl + 2: GOTO ILIM 


GOTO I LIM 

D1 = D1 + 2: GOTO ILIM 


ICURSVAL: ' *** FINDS THE VALUE OF THE LINE AT THE 
D1 = CSRLIN 

IF CSRLIN = 10 THEN GOSUB ICHANGE ' 
IF CSRLIN =12 THEN GOSUB ICHANGE ' 
IF CSRLIN = 14 THEN GOSUB ICHANGE • 
IF CSRLIN = 16 THEN GOSUB ICHANGE • 
IF CSRLIN = 18 THEN GOSUB ICHANGE ' 
IF CSRLIN = 20 THEN GOTO MENU • 
GOTO INSTRUMENT 


CURSOR LOCATION ************* 

ROI LOW 
ROI HIGH 
BKGND PERIOD 
DATA PERIOD 
CALIB . PERIOD 
RETURN MAIN MENU 


ICHANGE 

: 










FOR I = 

50 

TO 7 0 









NEXT I 

LOCATE Dl, 

I, 

0 

: PRINT 

II II 






IF 

IF 

IF 

IF 

IF 

Dl = 10 
Dl = 12 
Dl = 14 
Dl = 16 
Dl = 18 

THEN 

THEN 

THEN 

THEN 

THEN 

LOCATE 

LOCATE 

LOCATE 

LOCATE 

LOCATE 

Dl, 

Dl, 

Dl, 

Dl, 

Dl, 

60, 0: PRINT 
60, 0: PRINT 
58, 0: PRINT 
58, 0: PRINT 
58, 0: PRINT 

"(1 TO 
"(1 TO 
" ( 1 TO 
" ( 1 TO 
" ( 1 TO 

1024) " 
1024 ) " 
10, 000 
10,000 
10,000 

SECONDS) 

SECONDS) 

SECONDS) 

LOCATE Dl, 

49: INPUT 

SETTINGS 






IF Dl = 

10 

THEN 









END IF 

IF 

SETTINGS 

< 

1 

OR SETTINGS 

> 1024 THEN 

BEEP: GOTO ICHANGE 

IF Dl = 

12 

THEN 









END IF 

IF 

SETTINGS 

< 

1 

OR SETTINGS 

> 1024 THEN 

BEEP: GOTO ICHANGE 

IF Dl = 

14 

THEN 









END IF 

IF 

SETTINGS 

< 

1 

OR SETTINGS 

> 10000 THEN 

BEEP: 

GOTO ICHANGE 

IF Dl = 

16 

THEN 









END IF 

IF 

SETTINGS 

< 

1 

OR SETTINGS 

> 10000 THEN 

BEEP: 

GOTO ICHANGE 

IF Dl = 

18 

THEN 









END IF 

IF 

SETTINGS 

< 

1 

OR SETTINGS 

> 10000 THEN 

BEEP: 

GOTO ICHANGE 


instdata (Dl) = SETTINGS ' RESETS INSTRUMENT 

RETURN 


* * 


BACKGROUND : 7 ************************************************** AAVkJkAAA)|[ 

' *** BACKGROUND ROUTINES ************************************ 

BKGNDFILE$ = "BKGND. DAT" 



COSUB FLUSH 


“S£e !L 2 ?i. p ^?L , ':*r, BACKG ?? i ; N ?. Ro “TiN E s ...... 


LOCATE 11, 11; PRINT 
LOCATE 13, 11: PRINT ' 
LOCATE 15, 11: PRINT 1 

LOCATE 17, 11: PRINT " 
LOCATE 22, 14: PRINT 

D1 = 17 


BAKDLIM : 


[ J RECALL BACKGROUND DATA FROM DISK" 

[ ] SAVE BACKGROUND DATA ON DISK" 

[ ] ACQUIRE BACKGROUND DATA" 

[ ] RETURN TO MAIN MENU" 

USE <UP> OR <DOWN> ARROWS THEN <ENTER> TO SELECT." 


8" THEN D1 = D1 


MANAGES THE LIMITS OF TRAVEL OF THE CURSOR *** 

LOCATE > D1^ l3 EN l D1 0 = 13 ELSE IF 01 < 11 THEN 01 = 17 
D1 = CSRLIN ' ' ' 

" GOTO KSJT~ ELSE IF A * - “8” 

IF AS - "H" THE™™ i ? D I - I ™ T G0T0'BAKDLIM E G ° T ° BAKDLIM 

c^bak^m™™ 01 - - + - SSS bakdl™ 

D1 = D1 + 2 : GOTO BAKDLIM 
BAKCURSVAL: ' *** FIND o thp 

finds the value of the line at the cursor location ... 

01 = CSRLIN 

IF CSRT L ?m = }} THEN GOSUB BKGNDCALL 

IF CSRI?! = il THEN G0SUB BKGNDSAVE 

IF CSRT*TM = JS EN G0SUB BKGNDMEASURE 

if CSRLIN = 17 THEN RETURN 

GOTO BACKGROUND 

BKCNDCALL: - >*“:*;***“ 

' :** RECALL background data from disk 


RECALL BACKGROUND DATA 
SAVE BACKGROUND DATA 
MEASURE BACKGROUND 


GOSUB FLUSH 


bkgndfile$ 2 ^'.. bkgnd.dat** background Data recall ...... 

BKGNDCHK: ' ... CHECK FOR BACKGROUND DATA ... 

IF PNT = 1 THEN 

Se^i^n^t AEREA - ™ P “ GE " 

J? cl I Z. ° E 2 - -y THEN goto A s^^kLdfil S 
GOTO BKGNDCHK " """ ™ EN G0T ° BACKGROUND 

END IF 

IF r$ = "Y" OR 2 r $ X = ' " y " ^THEN °° ^ WANT T ° RETRIEVE A BACKGROUND DATA FILE (Y/N) 


GOTO SELBKGNDFIL 
ELSEIF r$ = "N" OR r$ = " n " THEN 
GOTO BACKGROUND 

END IF 

GOTO BKGNDCHK 

SELBKGNDFIL: ' *** SELECT BACKGROUND DATA FILE ON DISK *** 

GOSUB FLUSH 

LOCATE 12 2 19: P pR?NT^*R^vrS CKGR ° UND DATA RECALE ****" 

LOCATE 22, 16: PRINT "<ENTER> U FOR D DEFAm DEFAULT “ "» BKGNDFILES ; » ] . 

LOCATE 14, 48: INPUT ; A$ AULT FILE NAME OR TYPE IN OTHER NAME" 

IF A5 <> »» THEN BKGNDFILES = AS 
FOR B J G i! D 25 L TO 7 5 "“ THEN BEEP! GOTO SELBKGNDFIL 

NEXT J “ CATE 14 ' J ' 01 PRINT ” " ' CLEAR LINE 14 

FOR J = 15 TO 75 

NEXT J L ° CATE 12 ' J ' ° : PRINT " " ' CLEAR LINE 12 

CLOSE E U 2 ' 19 : PRINT " BACKGR °UND DATA FILE [ DEFAULT - " ; BKGNDFILES; " ] "• 

OPEN BKGNDFILES FOR APPEND AS #1 
XIMPUTFILE - LOF(l) 

CLOSE HI 

IF XIMPUTFILE <= 2 THEN 
KILL BKGNDFILES 

FOR P j i°2 A TO 78 ' 2< * ’' PRINT BKG NDFILE$; » NOT FOUND ON DISK" 

NEXT J “ ' J ' ° : PRINT ' ■ ' CLEARS COMMAND LINE 

LOCATE 22, 29: PRINT "HIT ANY KEY TO CONTINUE" 

LOCATE 5, 71: PRINT TIMES* 11 " , 

LOOP WHILE INKEYS = »" UPDATE REAL TIME 

GOTO BKGNDCALL 

END IF 

LOABKGNDFIL: ' ... LOADS BACKGROUND DATA FILE ... 

OPEN BKGNDFILES FOR INPUT AS #1 
FOR U I l 1 ! B £ G ” DDATE $> BKGNDTIMES, BKGNDTIME 
INST = (6 + (2*1)) 

NEXT I INPUT #1 ' instdata ( INST) 

FOR I = i TO 1024 

NEXT I INPUT * 1 * fa ackgroundarray& (I) 

instdata (14) ** BKGNDTIME 
CLOSE #1 


BKGND = 1 


' SETS BACKGROUND PRESENT FLAGS 


RETURN ' ********************************************************^^^^^^^ 
BKGNDSAVE: ' *** ***^** ************************************************* ^ 

; ;** SAVE BACKGROUND DATA TO DISK **************************** 

GOSUB FLUSH ********************************************************* 

IF BKGND =* 0 THEN 

WCATE L ?9 AT L-' “print "Sn T ANrKFj°TA A S^ R ° UND C 0 LLE C™ OR INSTALLED ... 
D0 iy, 24. PRINT HIT ANY KEY TO RETURN TO BACKGROUND MENU" 

LOCATE 5, 71: PRINT "••? TIMES ; / UPDATE RFAr ttmp 

LOOP WHILE INKEY$ = uPUAlt REAL TIME 

GOTO BACKGROUND 

END IF 

DSSK: ' *** DISK OUTPUT ROUTINE *** 

BEEP 

GOSUB FLUSH 

BKGNDFILID: ' *** FILe IDentif ication *** 

LOCATE 22 2 th BACKGR0U ND DATA SAVER ****•• 

Ti-ti v-™- sss oata <v». * 

JotS $ dssk D$ = " y " THEN G0T0 °wtdefault 


OWTDE FAULT: ' *** OUTPUT FILE DEFAULT NAME *** 


FOR A jV42 TO : ? 5 RINT " BACKGR0UND DATA FILE [ DEFAULT - ; BKGNDFILE$ ; •• ]"; 

LOCATE 14, J, 0: PRINT " " 

NEXT J 

“?a?e 14 : Hi Sgg Tlf ER> F0R DEFAULT ™ «ahe or type in other name .. 

IF A$ <> "i. THEN BKGNDFILES - AS 

FOR B J G = D 25 L ?o 75 " " THEN BEEP: G0T0 BKGNDFILID 

„_ v _ , LEGATE 14, J, 0 : PRINT " « 

NEXT J 

CL0SE E #1 2 ' 19 : PRINT " BACKGR0UND DA TA FILE [ DEFAULT = ; BKGNDFILES; " 

OPEN BKGNDFILES FOR APPEND AS #1 
XOUTPUTFILE - LOF(l) 

CLOSE U 

IF XOUTPUTFILE > 2 THEN BEEP: LOCATE 14, 29: PRINT BKGNDFILES; - ALREADY ON DISK 


OVWRT: ' *** OVERWRITE *** 


IF XOUTPUTFILE <= 2 THFN 

FOR J - 7 T0 T 75 N G0T0 SAVBKG NDFIL 

next J LOcate 22 , J, 0: PRINT « «• 

LOCATE 22 3? n* 

« F l? S B 'g^fi^s BS - ™en R MRITE (v ' n > BS 
end if NKGNDFILE$: goto savbkgndfil 

IF B$ = i'n” or n« „ i 

FOR 7 R - B ?c = n THEN 
r<JR j a 25 TO 75 

NEXT j L0CATE 14 / j, 0: PRINT '• H 
END IF G0T ° DSSK 

BEEP: GOTO OVWRT 


CLEAR COMMAND LINE 


CLEAR LINE 14 


*** 


SAVBKGNDFIL- ' *** 

Lj ' *** SAVEFILE 

NDDATE$ = DATE$: BKGNDIME$ • TIME$ 

open°bkgndfiles t FOR a ( 14 ) 

WBITE |i, BKGNDDl?L° U n^ T AS #1 

FOR I - 1 TO 6 ' BKg NDTIMES, BKGNDTIME 

WRITE*! ■ < 2 * I)) 

NEXT I TE #1, instdata ( INST) 

FOR 1=1 TO 1024 

next I WRITE bacJtgroundarrayfi (ij 
close #i 
RETURN ' 
bkgndmeasure: - 

, 4 ***********«**.»*.**...., 

' *** CONTROtc n» *************»,», 

; * ackground acquisition activities .... 

COSUB FLUSH 

BKSCREEN- • *** n ***************, 

-CATE 7 .‘ * BAC ~ ~ ... 

locate ii; « ; P ?lE :: ^ K =S“pERiSo Q ls rS - I0N M0DE ***” 

-CATE 22. 2, 0; PRINT^POMTIoirpIoBE 31 ™^ ^»»»»"' , ' SEC0NDS • , 

10 locate . „ . . * HIT <SPACE BARi T0 «» ~ 


LOCATE 5, 71, 0: print TIME$; „„ 


UPDATE REAL TIME 


A$ = INKEYS 

IF A$ = CHR$ ( 27 ) THEN GOTO MENU 
LOOP IF A$ = CHR$ (32) THEN GOTO BKOK 

BKOK: 

FOR I = 2 TO 78 

LOCATE 22, I, 0: PRINT " " 

NEXT I 


INTERRUPT 


' CLEARS COMMAND LINE 


1 

LOCATE 22, 18, 0: PRINT "MEASUREMENT IN PROGRESS, WAIT FOR COMPLETION" 
TIME% - instdata(14) 

GOSUB NORLAND 

' COME BACK FROM NORLAND 

MOVE DATA INTO BACKGROUNDARRAY 


FOR I = 2 TO 78 

LOCATE 22, I, 0: PRINT " " 

NEXT I 


7 CLEARS COMMAND LINE 


LOCATE 22, 17, 0: PRINT "DATA TRANSFER IN PROGRESS, WAIT FOR COMPLETION" 
FOR x% = 1 TO 1024 

NEXT x% backgroundarra y & ( x %) - chanreadS (x% , &HF300, &HD000) 


BEEP 

FOR I - 2 TO 78 

NEXT I L0CATE 22 ' J > 0: PRINT " » 


CLEARS COMMAND LINE 


LOCATE 22, 16: PRINT "MEASUREMENT OVER, HIT ANY KEY TO RETURN TO MENU" 
DO 

LOOP WHILE C INKEVS PRINT TIMES! "" ' UPDATE REAL TIME 


BKGND = 1 


BACKGROUND DATA FLAG 


GOTO MENU 

t *** CONTROLS DATA ACQUISITION ACTIVITIES ************************* 




' GOSUB FLUSH 
IF BKGND = 0 THEN 

LOCATE & Hi 

LOCATE 5, 71, 0: PRINT ""; TIMES- "•• 

LOOP WHILE INKEYS = "" 

GOTO MENU 
END IF 

IF CAL = 0 THEN 

?0 S RE?u^ F ?fJ“™Lu D “ A -" 

Trtrin ,¥J CATE 5 ' 71 ' 0: PRINT TIMES- •"» 

LOOP WHILE INKEYS = " " ** 

GOTO MENU 
END IF 

MEASURESCREEN; ' ... MEASUREMENT SCREEN DISPLAY ... 

““?f ?i 26 2 i.T B rL"::L^ TA ACQUISITION MODE .... 

LOCATE 22, 3, 0: PRINT u POSITION T pROBE OD THEN'HIT nStdata ( 16) ’’ " SECOf, DS" 

D0 R0BE ' THEN HIT <SPACE BAR> TO BEGIN MEASUREMEN 


LOCATE 5, 71, o: PRINT TIMES* / nonxmo 

nt *' UPDATE REAL TIME 

' INTERRUPT 


A$ = INKEYS 

Jp A f CHR$(27) THEN GOTO MENU 
LOOP IF AS ~ CHR$ ( 32 ) THEN GOTO MSOK 


MSOK: 


LOCATE 

LOCATE 

LOCATE 

16, 52: 
18, 52: 
14, 50: 

PRINT " 
PRINT " 
PRINT " 

II 

II 

FOR I * 

2 TO 78 



NEXT I 

LOCATE 

22, I, o: 

PRINT 


CLEARS COMMAND LINE 


LOCATE 22, 18, 0: PRINT "MEASUREMENT IN PROGRESS, WAIT 

TIME% = instdata (16) ' DATA ACQUISITION PERIOD 

GOSUB NORLAND 


FOR COMPLETION" 


COME BACK FROM NORLAND 
MOVE DATA INTO DATAARRAY 


FOR I = 2 TO 78 

LOCATE 22, 

NEXT I 

LOCATE 22, 17, 0: 


1/ 0: PRINT " " 


' CLEARS COMMAND 


PRINT "DATA TRANSFER IN PROGRESS, WAIT FOR 


LINE 

COMPLETION" 



FOR x% 


= 1 TO 1024 

dataarray& (x%) = chanread& (x% , &HF300, &HD000) 
NEXT x% 1 


NORMALIZE BACKGROUND COUNTS PER CHANNEL TO DATA PERIOD 
SUBTRACT NORMALIZED BACKGROUND FROM DATA 
PUT PROCESSED DATA IN PROCESSEDARRAY 
PROCESSED ARRAY IS DEMENS IONLESS IN TIME 


FOR I = 2 TO 78 

NEXT I L ° CATE 22 ' T ' ° : PRINT " " ' CLEARS COMMAND LINE 


LOCATE 22, 18, 0:- PRINT "CALCULATION IN PROGRESS, WAIT FOR COMPLETION" 
FOR x% = 1 TO 1024 


processedarray ! (x%) « 
IF processedarray! (x%) 


(dataarrayfi (x%) - (backgroundarrayfi (x%) 
< 0 THEN processedarray! (x%) = o 


* (instdat 


NEXT x% 


PERFORM AVERAGE CHANNEL CALCULATION 
weightedcounts& = o 

COUNTSTOT& = 0 ' ZERO OUT THE TOTAL 

FOR I = instdata (10) TO instdata(12) ' ROI LOW TO ROI HIGH 
COUNTSTOT& = COUNTSTOT& + processedarray! (I) 

weightedcounts& = weightedcountsfi + (processedarray! (I) * i ) • CHANNEL FO 

NEXT I 


AVECHAN! - 
Zeff = (-B 
LOCATE 14, 
LOCATE 16, 
LOCATE 18, 
'LOCATE 19, 
BEEP 
BEEP 

' LOCATE 20, 
' LOCATE 13, 
' LOCATE 15, 


weightedcounts& 
+ SQR(B * B 
23: PRINT " 

23: PRINT " 

23: PRINT " 

23: PRINT 


23 

23 

23 


/ COUNTSTOT& 

4 * A * (C - AVECHAN!))) 
TOTAL COUNTS IN ROI = » « 
AVERAGE CHANNEL NUMBER * 
EFFECTIVE ATOMIC NUMBER 
A ; B ; C 


/ (2 * A) 
COUNTSTOT& 

■ "; USING "####•# 1 

■ USING 


PRINT 

PRINT 

PRINT 


"time/tyme ="; instdata (16) / instdata (14) 
"Lowest channel# =" ; instdata (10) 

"Highest channel# =»; instdata (12) 


; AVECHAN! 

Zeff 


FOR I 


2 TO 78 

LOCATE 22, I, 0: PRINT " " 


/ 


CLEARS COMMAND LINE 


ne^t I 


GOTO MEASURESCREEN 


get ready for next measurement 


R F'Ttidkt / A UKcJMENT 

elush: . 

; :i”i” :i ” henus md " pdates ° ate «"«««**■ 

Bin k r-» ««*, .... 

NEXT I L ° CATE J » 3 / 0: PRINT '• 

FOR J = 2 TO 78 

NEXT J L0CATE 22 > J/ 0: Print ii ii 


’ UPDATE REAL TIME 
CLEARS MENU AREA OF SCREEN 


CLEARS COMMAND LINE 

cZr . ZZZ!ZZ‘*‘******‘**’‘‘‘‘'‘*‘’‘*****‘*** 

»„:;;**************************- 

18< 0: print data -" 

LOOP WHILEINKEYs'.°.m| PRINT "" ; TI «E$ ; »» 

GOTO MENU 

SND IF 


,*** MEASUREMENT SCREEN DISPLAY *** 


•OCATE 12, 2 ii ; P pRiJ T "*** CALIBRATION MODE ***>• 

•OCATE 14, II; PRINT •• f ™L£ EFAULT CALIBRATION" 

oca?e e it : ii; Si5 : d i^^ LIBRATION constants- 

D1 . 12 <UP > ° R <“WN> ARROWS THEN <ENTER TO SELECT » 

SSHvl »-V<- 

D1 = CSRLIN ' ' 13 

IF A$ r - KE 2 ” ' T HEN A Dl"= C DI 1 GOTO Si"®*' 1 ELSE IF « = "8" THE! 


IF A$ N i A »H"~THEN H Di A$ n 7 RIGHT $( A $'. X ) ELSE GOTO MLIM1 

TP £o - ™ EK D1 - D1 - 2: GOTO MLIM1 

GOTO MLIM1 01 = 01 + 2 : G0T0 MLIM 1 

D1 = D1 + 2: GOTO MLIM1 

MCURSVAL1: ' **‘ l FINDS THE VALUE OF THE LINE AT THE CURSOR LOCATION ............ 

IF gSS': “ THEN S NEHCON T ° SE DE ™ E * CALIBRATION 

s^sss - 16 ™ en ™ i'Ssr* 5 *■ b ' »• 

RETURN 

DEFAULT: 

GOSUB FLUSH 

LOCATE 13,' 31: PRINT "A^-q^ 06" rBRATI0N CONSTANT VALUES" 

LOCATE 15, 31: PRINT »B » 14 ! 072" 

LOCATE 17, 31: PRINT "C =» 27.189" 

CAL = l 22, 18 PRINT "HIT <ESC> TO RETURN TO MAIN MENU." 

DO 

A$ = INKEY$ 

IF A$ = CHR$(27) THEN GOTO MENU 

TnnD CATE 5 ' 71 ' 0: pr INT TIMES- »" 

LOOP UNTIL A$ = CHR$ (27 ) * 

RETURN ; 

NEWCON: 

GOSUB FLUSH 
LOCATE 12, 11: 

LOCATE 14, 11: 

LOCATE 16, 11: 

CAL = 1 

LOCATE 22, 18: 

DO 

A$ = INKEYS 

loCATE = 5 CH ?f <2 n ) ™ EN G0T ° MENU 
t^^ CATE 5 ' 71 ' 0: PRINT TIMES- "" 

LOOP UNTIL A$ = CHR$ ( 27 ) 

RETURN V ; 


INPUT " A 
INPUT " B 
INPUT " C 


? 

« 


A 

B 

C 


PRINT “HIT <ESC> TO RETURN TO MAIN MENU.” 


RECAL: 

GOSUB FLUSH 

LOCATE 10, 22: PRINT "MEASUREMENT PERTon to ii 

LOCATE 22, 3, 0: PRINT "HIT <ESC> TO USp ™»firi nStdata (18) ; "SECONDS" 

DO <ESC> T0 US E DEFAULT CALIBRATION. HIT ANY OTHER KEY 

A$ = INKEYS 

mA!p S e CliR ^ 27 ' THEN GOTO MENU 
t EOCATE 5, 71, 0: PRINT TIMES' "" 

LOOP WHILE AS ■ "" 

FOR J = 2 TO 78 

LOCATE 22, J, 0: PRINT " » 

NEXT J 

LOCATE 13, 3, 0: INPUT " 

FOR M = 1 TO N 


NUMBER OF CALIBRATION POINTS = ? (MUST BE >2) 


• LOCATE 20, 8 
LOCATE 17, 3 
LOCATE 22, 3, 0: 

DO 


n! THE REMAI NING NUMBER OF CALIBRATION 

0. INPUT " EFFECTIVE 2 OF TARGET IS *> » 

RINT "POSITION PROBE, THEN HIT <SPACE BAR> TO 


POINTS IS " ; N 
Z(M) 

BEGIN MEASUREMEN 


LOCATE 5, 71, 0: PRINT ; TIME$ ; " " 
A$ = INKEY$ 

IF A$ = CHR$ ( 27 ) THEN GOTO MENU 
.... IF " CHR$ (32) THEN GOTO CALOK 


' UPDATE REAL TIME 


INTERRUPT 


CALOK: 


FOR I = 2 TO 78 

LOCATE 22, 
LOCATE 17, 

NEXT I 


I, 0: PRINT " •' 
I, 0: PRINT " » 


CLEARS COMMAND LINE 


LOCATE 22, 18, 0: PRINT " 
TIME% = instdata(l8) 

GOSUB NORLAND 


MEASUREMENT IN PROGRESS, WAIT FOR COMPLETION 
' DATA ACQUISITION PERIOD 


II 


COME BACK FROM NORLAND 
' M0VR data INTO DATAARRAY 


FOR I = 2 TO 78 

NEXT I L ° CATE 22 ' I ' ° : PRINT " " ' CLEARS COMMAND LINE 

LOCATE 22, 17, 0: PRINT "DATA TRANSFER IN PROGRESS, WAIT FOR COMPLETION" 
FOR x% = 1 TO 1024 

NEXT x% dataarray& ( x% ) = chanreadfi (x% , &HF300, &HD000) 


' gS E NO B ^L?^rB“KGRO U( f D R F^of T L A TO ““ PERI ° D 

FOR I = 2 TO 78 

LOCATE 22, I, 0: PRINT " " 

NEXT I 1 


' CLEARS COMMAND LINE 


LOCATE 22, 18, 0: PRINT "CALCULATION IN PROGRESS, WAIT FOR COMPLETION" 
FOR x% =1 TO 1024 


processedarray ! (x%) = 
IF processedarray! (x%) 


(dataarrayfi (x%) - (backgroundarray& 
< 0 THEN processedarray! (x%) = o 


(x%) 


* (instdat 



NEXT X% 


' PERFORM AVERAGE CHANNEL CALCULATION 
weightedcountsSi = 0 

COUNTSTOT& = 0 ' ZERO OUT THE TOTAL 

FOR I = instdata ( 10) TO instdata(12) ' ROI LOW TO &>I HIGH 
COUNTSTOT& - COUNTSTOT& + processedarray ! (I) 

weightedcounts& = weightedcountsSi + (processedarray! (I) * i) • CHANNEL FO 


NEXT I 

AVECHAN! = weightedcountsSi / COUNTSTOT& 
E(M) - AVECHAN! 

BEEP 

FOR I = 2 TO 78 

LOCATE 22, I, 0: PRINT " " 

LOCATE 17, 3, 0: PRINT " M 
NEXT I 


NEXT M 

' Least-Squares Calculation Routine . 

program to calculate the coefficients 
'fit of a second-order polynomial: 

' <chj|> = A*Z A 2 + B*Z + C 
'to calibration data points. 


for a least-squares 


'Least-Squares Calculation 


SI - 0: S2 = 0: S3 = 0: S4 = 0: S5 


FOR I = 1 TO N 

51 - SI + Z(I) 

52 = S2 + Z(I) * 

53 = S3 + Z(I) A 

54 - S4 + Z(I) A 

55 - S5 + E (I) 

56 =» S6 + Z(I) * 

57 - S7 + Z(I) * 
NEXT I 


Z(I) 

3 

4 

E (I) 

Z(I) * E ( I ) 


0: S6 


0: S7 = 0 


\ 


^ * S4 — S3 * S7\ 

B2 = S5 * fs? it ci » 

B3 = S2 * (sl * S7 " ^ * S4) 
(bl * S7 - S2 * S6) 

B = (B1 + B2 + B3) / D 
A (A1 +. A2 + A3 ) / d 


Cl = S5 * 
C2 = si * 
C3 *■ S2 * 


(S2 * S4 - S3 * S3) 

S3 * S7 - S4 * S6) 

(S6 * S3 - S7 * 82) 


C = (Cl + C2 + C3) / D 

FOR I = 2 TO 78 

LOCATE 20 ' i' ° ! PRINT " " 

NEXT I /I/O: PRINT 11 m 

“cate iSJfpggV 

“il £ S 3 * :«?■*«« co EPFIC1ENTS ,. 

sS : SJ S "c B : ::l 

LOCATE 12 * 50 : p^ NT (TRUE) •' 

FOR I a i'to n RINT " 2 (CALCULATED) " 

LOCATE 13 + 40 •* PRINT 2(1) 

NEXT I ' *' 50: PRI W (-B + SQR(B * B - 4 * a * (c 


£ (I)))J / (: 


CAL = l 

Jp” 22, 3# 0 : PRINT ..CALIBRATION IS COMPLETE PR ESS <ESC 

A$ * INKEY$ PRESS <ESC> TO RETURN 

LOOP^ ch R$(27) THEN GOTO MENU 


RETURN 

NORLAND: ' ***********„„ 

' *** NORLAND MOi ***************************, 

r ° umnes 

' p HA Initialization: ********“*“*****‘*************‘****....*.,,, 

' offset^T*!^ *^ r ?HA n“e? r ADD et DI f R r fj! 00h ! IRQ7 > D °°°°h 

gain! . 102 ; S6C ” ebasa (LIVEti “>.'^ slL: 20 * Uld - 



CALL setoffset (0, gain%) 
CALL setlld(20) 

CALL setuld(1023) 

OUT &HF307 , &HD5 


set 0 offset and lk gain 
' S ET LLD TO CHANNEL 20 
' SET ULD TO CHANNEL 1023 

now write inode register 


; * mp mode _ bi t 7 =o, PHA mode - bit 6 set 
, l°«er memory bank - bit 5 =o, ADD - bi?J'set 
1st lk memory, bits 3-o are 0101 4 set 

FOR AMP MODE USE &H55 INSTEAD OF &HD5 


OUT &HF304 , &H8C 


' 100ms base rate * lo 


““ SE; s s i£So°; SHD00 °> ; f££*i££r* ac 

) acquire TIME% sec realtime 

' for REAL TIME USE CALL acquire (TIMEt, SHF0, RHF300) 

“ 

FUNCTION chanread (A%, b% c%i 
OUT B% , &H41 ' 

DEF SEG = c% 

tempo = PEEK (A% * 4 ) 
tempi - PEEK(A% * 4 + 1 ) 
temp2 = PEEK (A% * 4+21 
chanread = temp 2 . 65536 + tempi . 256 + tempo 

END FUNCTION 

SUB setlld (lld%) 

SHARED gain% 

: Registers^^nd 1 " a?e affec^ed^H leVSl disc ^minator 

are affected and must be rewritten after this routine 

' ValUS WrUten tQ re ^isters depends on conversion gain 
lld% = Hd% * 4096! / gain% 


upperbyte% = FIX(lld% / 256 ) 
lowbyte% - FIX(lld% MOD &HFF) 
OUT &HF307 , upperbyte% } 

OUT &HF300, &H80 
OUT &HF301, lowbyte% 

OUT &HF300, &H4 0 
END SUB 

SUB setoffset (offset%, gain%) 


' separate upper byte 
' separate lower byte 

' enahle P ?fa 4 bi ? s of lld in Register? 

' write to^ld*;! 113 * 6 ? and dlscrim inator buffs 
write to lld register. _ The upper 4 

' now disable discriminate? ??gis???s ar * alS 


9 - 


V 

, > 

' This routine will set the proper offset and will adjust 
factory offset for different gains 
' The conversion gain is also set by this routine 

DEFLNG V 


' First get factory offset 
IF gain% = 8192 THEN 


' A note about the hexadecimal numbers: 

: chop o« t the e feldi^°n?° f ^ 2 f ? r a l0ng inte 9*r. QuickBasic will 
, . ° p Off the leading 0's when i move to the next line. When it 

interprets the number, it is sign extended, leaving ffffff42 which 

: "6«: fi So a L have ?. al ^ er 2 at f ves: ^ S h 

decimal, 65346, and specify in decimal or 2. Do somethina like 

, i[Je S Lin 1 S 00ff4 ^ AND &H2000ff42 which will be interpreted as 0000ff42 

' the fi?st me?Li S w??rh re ? dable f rom a P ro< ? ra mming point of view, but 
the first method will be faster in execution speed. 

VALUES = SH1000FEE9 AND SH2000FFFF 

prcpTr '.factory offfset 3fe9 + gain bits 11 
ELSEIF gam% = 4096 THEN 

VALUES - SH1000BF65 AND SH2000FFFF 

ftcptp ' . fact ory offset 3f65 + gain bits 10 
ELSEIF gain% = 2048 THEN 

VALUES - SH7FA4 

' factory offset 3fa4 + gain bits 01 

ELSEIF gain* n ^l55r?HEN n ° Si9n bit ' the 7tba is rea<3 “meetly 
VALUES = SH3FC3 

END IF # factory offset 3fc3 + gain bits 00 


' N°w subtract desired offset 
VALUES = VALUES - offset% 

' now write the registers 

lowvalue% = VALUES AND SHFF 
hivalue% = FIX (VALUES / 256) 

OUT SHF302 , lowvalue% 

OUT SHF303, hivalue% 

OUT SHF305 , 0 

END SUB 


' get lower byte 
' separate upper byte 

' output offset (0 + factory offset for 8k) 
' upper offset 3FH + gain bits 2K 
load the offset into the hardware 


SUB setuld (uld%) 
SHARED gain% 


/ 

/ 

/ 

t 


This routine will set up the upper 
Registers 0 and 7 are affected and 
value written to registers depends 


level discriminator 

must be rewritten after this routine 

on conversion gain 


ul«% = uld% * 4096! / gain% 

upperbyte% = FIX(uld% / 256) 
lowbyte% =■ FIX(uld% AND &HFF) 
OUT &HF307, upperbyte% 

OUT 6HF300, &H0 
OUT &HF301, lowbyte% 

OUT &HF300 , &H40 




* 


separate upper byte 

separate lower byte 

put upper 4 bits of uld in Register? 

enable uld register and discriminator buffe 

write to uld register. The upper 4 

' bits from reg7 are als 
now disable discriminator registers 


END SUB 
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Abstract 

A compact Instrument capable of determining the 
etfecuve atomic number (Z^) of a target has been 
developed and characterized. It is a non-destructive sensor 
which utilizes a cadmium telluride (CdTe) solid state 
dealer ana a »Sr/»Y aourca in a naakacaiinr geomeln, 
The spectrum of the backscanered beta panicles depends 
on the target atomic number. Measurements show that by 
using the mean energy of the backscatter spectrum, the 
average effective atomic number can be determined to 
better than ± l Z. 


Introduction 

In many applications which use ionizing radiation as a 
probe, the size, weight and power requirements of the 
nuclear sensor are important considerations, particularly 
for portable Instruments and instruments designed for use 
In space. With this in mind we undertook the development 
of a radiation-backscatter sensor system which is small 
lightweight and uses little power. The design of the system 
can be adapted for a variety of non-destructive analysis 
applications by selecting the type and energy of the 
radiation source used. The system described here was 
designed to determine the average Z of a target by 
measuring beta backscatter. 7 

Beta backscatter is largely caused by coulomb scattering 
from the target nucleus, (1,2]. As a result, the spectrum of 
backseat tered betas depends on the target atomic number. 
Tnts was demonstrated in an early experiment by Brownell 
[3] who showed that the average energy of the 
backscanered spectrum increases monotomically with 
target Z. We have exploited this relationship to determine 
inc of target materials. 

The system developed includes a sensor head (Figures 1 
and 2), electronics and software. The sensor head contains 
both the radiation detector and a radioactive isotope. The 
sensor chosen is a CdTe particle detector we manufacture. 
The source is a 2 mCi *>Sr/*Y ceramic pclleu This 
ceramic source is rugged and safe, preventing potential 
contamination. The maximum energy of the betas is 2.3 
MeV, 

The source is embedded into a threaded plastic piece 
which screws into the sensor head, allowing for adjustment 


of the beam collimation. The beam collimation if 
adjustable from a beam spread of 15* to 50*. The field oi 
view of the detector is larger than the area illuminated by 
the radiation at any collimation, which allows the beam 
collimation to define the target area of analysis. 


*OU«C* COLLUhTO* 



Figure 1. Diagram of sensor head for beta backscatter 
measurements. 



Figure 2. Photograph of sensor head for beta baclcscaiii 
measurements. 
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Because of ihe close proximity of ihc source and 
deieclor, background noise from the source is a concern. 
Several steps were. taken to minimize the background while 
maintaining the desired small size. When using beta 
emitting sources, background noise is due to both the 
direct beta particles themselves and to the brcmslrahlung 
radiation generated as the particles pass through matter. 

The elimination of the noise due to the direct detection 
of the beta particles simply required that the minimum 
wll thicknesses be equal to or greater than the maximum 
range of the betas. The wall thicknesses will vary 
depending on the material used. This, however, is not the 
most important source of noise. 

A more difficult problem stems from the fact that beta 
panicles generate bremsslrahlung radiation as they p ass 
through matter. The magnitude of the bremsslrahlung flux 
produced is a linear function of Z; minimizing Z minimizes 
ihe Ilux of bremstrahlung radiation. Thus, the sensor head 
is made from a low Z plastic. There is still a significant 
amount of background from the brcmslrahlung generated 
bp the plastic. Reducing this to acceptable levels required 
ihat about 4 mm of lead be added between the sensor and 
source. 

A final consideration in the choice of material used to 
make the sensor head is that many plastics degrade under 
mnstanl exposure to radiation. To address this, Dclrin™ 
plastic was chosen because It is relatively insensitive to 
radiation. 

Electronics 

The electronic circuitry for the system consists of an 
RMD, Inc., PCD preamplifier and an RMD, Inc., KS2 PC- 
board linear amplifier with bias supply which is inserted 
itto the PC The output of the linear amplifier is sent to a 
Norland 5000™ PC-board MCA. The linear amp, 
detector bias supply, and MCA are all installed in a' 
portable personal computer. The detector signal is 
toupled to the charge sensitive preamplifier via six feet of 
tmall diameter, low noise shielded cable. 

The computer controls the calibration, data acquisition, 
ind data analysis. A series of menus leads the user through 
calibration, collection of a reference background and 
collection of data. The MCA stores the data in RAM. 
Background is subtracted from the spectrum, then the data 
b put through a simple weighted mean algorithm. This 
produces a mean channel number (each channel represents 
is energy bin). The calibration routine does this for 
ttveral known samples and calculates the constants for a 
iklrd order polynomial fit. 

To analyze a target material, the sample is scanned and 
the mean channel number from the unknown target 
compared to the calibration curve. Based on the result a 


Zgfl value is assigned to it. This is displayed to the user 
along with an uncertainly value based on the total 
integrated counts from the measured spectrum. In the 
final design, a look-up table will be added, to identify likely 
materials which fit these parameters. 

Experimental Set-up 

In operation, the face of the sensor head is positioned 
so that the collimated beam irradiates the target area of 
interest. The dose proximity of the detector to the source 
gives a backscaucr angle of approximately 180°, as long as 
the target distance is greater than the minimum defined by 
the system. The cone of the collimated beam and the 
sensor head to target distance determine the area of 
interrogation and the flux of backseat lered betas viewed by 
the detector. 

In order to get reproducible results, there are several 
factors which need to be considered before making 
measurements. The radiation beam takes the form of a 
cone, which irradiates an increasing area with increasing 
sensor-to-target distance. The target must be large enough 
to fill the entire field of Interrogation at the distance and 
collimation being used. The targets may have uneven or 
irregular shapes, but must be of uniform composition. 
Finally, the targets must effectively be ’infinitely thick* to 
the radiation. According to measurements made by 
Sharma and Singh (4) and confirmed by measurements 
made in this study, infinitely thick is approximately equal 
to 20% or the maximum range of the maximum energy of 
the beta emissions being used. For the ^t/^Y source 
being used (maximum energy - 2.3 MeV), this translates 
to less than I mm for targets with densities greater than 1.5 
g/cm- . All targets used in this study were infinitely thick. 

Another factor which must be taken into account is the 
effect of air on the beta particles. Although the range of 
the 2.3 MeV betas is more than 7 meters in air, a 
significant fraction of the lower energy betas can 
experience coulomb scattering and energy loss through 
ionization in distances as short as a few centimeters. This 
changes the backscattcrcd spectrum as a function of 
distance. Additionally, scattering of the lower energy betas 
in air causes an expansion of the beam cone. Targets which 
otherwise would be large enough to fill the field of the 
beam may be too small as a result of the scattering of low 
energy betas in air. When a target is too small, the low 
energy betas which have air scattered to the outer edges of 
the beam miss the target and are not backscattcrcd. This 
loss of low energy betas from the backsca tiered signal shifts 
the mean energy upwards. 

Experimental Results 

The CdTe detector was operated with the negative 
contact irradiated. Since the beta particles do not 
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mean energy was determined and plotted as a function of 
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less than about 3.5 cm, angular effects cause the 
backseat, cr angle to deviate sufficiently from 180“ to alter 
the backscattcr spectrum. Therefore, 3.5 cm was used a 
the minimum operating distance of the device. 
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atomic number of a target material 
can energy of backscattered beta 


Figure 4. The sensitivity of the beta backscattcr device to the Zetr 
Is shown to be better than x 1 Z. ** a 


Figure 5. The effect of varying the detector bias on sensitivity tc 
Atf. 24 volts was chosen as the operating bias. 

The effective 
determines the 

particles. Measurements were made to determine t k 
sensitivity of the sensor head to changes in Z- n . Severi 

“1 Wdl 35 bina,y compounds wJre tested 

Re^Tii. C uncertain ‘y °' Zcff «o be about xl Z 

Repeated measurements indicate that the uncertainty ia 
Uie aknlaied Z^of the target is slightly larger than N ** 
*L ,hc louI integrated counts. This holds for all N 
,‘ he . Z e' r unccrla inty reaches near 1. Increasing the 
Strain n o‘ improve the accuracy significantly. 
Extrapolation of the data gives an ultimate accuracy of 


TTic sensor system performance was characterized by 
measuring the following parameters: ^ 


1) mean energy as a function of sensor-to-tareet 

distance 6 

2) mean energy as a function of 

3) mean energy as a function of density 

4) the effect of varying detector bias voltage. 


Measurements were made to determine the effect 
bfas on the mean energy calculation being used 
determine Z^. CdTe detectors can operate over a wk 
range of bias voltages. For the 1 mm device chosen f< 
these measurements, the bias could be varied from about 
volts 10 more than 60 volts. Higher bias imptovi 
sensitivity to higher energy radiations but increases th 

™ lsc : ,7*® rcsuIts («« F<Cure 5) indicate that above 
threshold of about 20 volts, there is no net effect of bias o, 






ik determination of Z^. Therefore, 24 volts was chosen 
a ihc operating bias since this is above the threshold and 
ikis voltage is readily available in most applications. 

The influence of density on the backscattcr signal was 
>bo examined. The oxides of beryllium and hydrogen shed 
«ful light on this question, since they have Identical Z- n , 
but the beryllium compound is 2.9 times more dense than 
ike water. The mean energy of the baclcscattered betas is 
virtually identical for both targets. Moreover, the trio of 
graphite, silicon, and potassium chloride arc of very similar 
densities; 2.2, 233, and 1.98, respectively. However, the 
ncan channel numbers of their spectra, 34.2, 59.4, and 
J6.2, mirror the progression in their Z^: 6, 14, 18.3. 

These observations suggest the conclusion that density 
toes not influence the mean energy of the backscattered 
betas. Additionally, the smoothness of the curve in 
Figure 4 supports this conclusion; since the densities of 
materials with increasing Z^ fluctuate dramatically. 

Conclusions 

A compact backscatter analysis system was developed 
Mscd on a CdTe detector and an isotopic beta source.' The 
system includes data acquisition and analysis software 
which runs on a personal computer. Testing indicates that 
ii was possible to distinguish the average atomic number of 
i material to ±1 Z regardless or the density of the 
material. Further improvements promise to give better 
than ±0.5 Z sensitivity. The system was designed to be 
flexible enough to be applied to a variety of applications by 
changing the data analysis algorithms and the source tvne 
md energy. 
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